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How to reconcile reversibility of microscopic laws with 
the arrow of time prescribed by the second law of thermo¬ 
dynamics? Ever since the pioneering work of Boltzmann, 
the question has stirred debate and controversy. With¬ 
out solving the core issue, Onsager realised that micro¬ 
reversibility has implications for macroscopic irreversible 
laws [H, Q. In particular, in the regime of linear response 
around equilibrium, micro-reversibility implies the sym¬ 
metry of properly defined Onsager coefficients. Building 
on preceding case studies , Brandner et al. Q devel¬ 
oped in a recent, remarkable paper the stochastic ther¬ 
modynamic formalism for periodically driven systems. In 
particular, they identify explicitly the Onsager coefficient 
describing the regime of linear response averaged over one 
period, and show, amongst other, that they obey a gen¬ 
eralized symmetry principle = Lji, where the tilde 
refers to the same system but driven by the time-reversed 
periodic signal. Under time-symmetric driving, the usual 
Onsager symmetry = Lji is recovered. They next dis¬ 
cuss the issue of operating without entropy production, 
hence reaching reversible efficiency. As was pointed out 
earlier for systems subjected to a magnetic field 0 , it 
appears that one could reach such an efficiency while op¬ 
erating at finite power whenever the Onsager matrix is 
not symmetric. However, a more detailed case study for 
a magnetic system 1T| and one in the Brandner et al. pa¬ 
per for a system subject to time-asymmetric driving show 
that “additional” restrictions appear, which prevent this 
from occurring, see also I2l4l6l|. It is speculated that 
these extra conditions lie outside the realm of the sec¬ 
ond law. In the present letter we show that this is not 
the case, provided the full impact of the second law in 
periodically modulated systems is assessed. More pre¬ 
cisely, one needs to take into account the symmetries im¬ 
plied by the second law on all the Onsager coefficients 
that appear, by considering the decomposition of the pe¬ 
riodic perturbation in its Fourier components. While, 
under time-asymmetric periodic driving, asymmetric On¬ 
sager matrices may appear in a coarse grained descrip¬ 
tion, the fine structure of the Onsager matrix implies 
that all fluxes, including the reduced fluxes, converge to 
zero in the limit of zero dissipation. The implication for 
any Onsager matrix, reduced or not, is that its determi¬ 


nant has to vanish in the limit of zero dissipation. In 
the particular case of only two fluxes, the corresponding 
Onsager matrix becomes symmetric. 

We start by reviewing the aforementioned thermody¬ 
namic puzzle. Consider a system that is brought out 
of equilibrium by the application of two thermodynamic 
forces Fi and F 2 , for example a temperature gradient 
and a chemical gradient. As a result, the system is no 
longer at equilibrium and corresponding fluxes Ji and J 2 
appear, for example a heat and particle flux. The com¬ 
monly considered situation is that of a non-equilibrium 
steady state, so that all these quantities are time inde¬ 
pendent. For proper choices of the fluxes and forces, this 
non-equilibrium steady state is characterized by a steady 
entropy production rate given by: 

S = FiJi + F2J2 > 0 . ( 1 ) 


The second law specifies that this quantity cannot be 
negative. Of special interest is the case in which an en¬ 
tropy decreasing loading process FiJi <0 (for example a 
particle flux uphill a chemical gradient), is made possible 
by an entropy producing driving process F 2 J 2 > 0. The 
second law becomes a statement about the efficiency ry of 
this transformation: 


11 = - 


FiJi 

F2J2 


< 1 . 


( 2 ) 


The maximum efficiency I is obviously reached for a re¬ 
versible process, i.e. a non-dissipative process with zero 
entropy production. For zero forces, the system is at 
equilibrium and the fluxes vanish. Hence a linear rela¬ 
tion is expected between fluxes and forces when the latter 
are sufficiently small: 


Ji = LiiFi + L12F2 

J2 = L21F1 + L22T2. ( 3 ) 

In this region of linear response, the entropy production 
reduces to the quadratic expression: 

S = L^iFf + (Li 2 + L21) F1F2 + L22FI ( 4 ) 

The second law, requiring the non-negativity of this 
quadratic form, gives rise to the following conditions: 

Lii > 0 , L22 > 0 AL11L22 > {Li 2 + T 2 i)^. ( 5 ) 
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Zero entropy production is trivially reached for Fi = 
F 2 = 0. The question of interest is whether this sit¬ 
uation, and hence efficiency 1, can also be realized for 
non-vanishing values of the forces. Since the quadratic 
form Eq. (|3]) is non-negative, this can only happen when 
its zero value is degenerate, i.e., when the discriminant 
is zero: 

4L11L22 = (Ti 2 + ^21)^- (6) 

Under this condition, the rate of entropy production can 
be rewritten as S = (^/LnFi ± \/T^F 2 ) , and thus van¬ 
ishes along the line F 1 /F 2 = L 22 /Ln. The cor¬ 

responding efficiency is equal to one. Since zero en¬ 
tropy production corresponds to equilibrium, a reason¬ 
able guess is that, under this condition, the fluxes should 
also vanish. With the above relation between the forces, 
one however finds that the fluxes: 

( Li 2 - L21) {L21 - L12) 

= - 2 - 2 - i"!, [7) 

only vanish for a symmetric Onsager matrix L 12 = ^ 21 - 
Insisting upon this reasonable line of thought, one con¬ 
cludes that, in the limit of zero entropy production, the 
2x2 Onsager matrix has somehow to become symmet¬ 
ric. Another related and revealing observation is to real¬ 
ize from the linear relation Eq. o, that vanishing fluxes 
can only appear for nonzero forces if the determinant of 
the Onsager matrix is zero, L 11 L 22 = T 12 T 21 . This, com¬ 
bined with a zero discriminant, cf. Eq. again implies 
Li 2 = ^ 21 - A zero determinant means that the fluxes are 
proportional to each other, a property which has been 
called strong coupling. The above discussion prompts us 
to ask the following questions. 1) How do we reconcile 
the above conclusion with the fact that 2x2 Onsager 
matrices can be asymmetric? 2) Is the extra condition, 
rendering the Onsager matrix symmetric in the zero dis¬ 
sipation limit, extraneous to the second law since it does 
not follow from Eq. ([5|)? 3) What about larger Onsager 
matrices? We are able to answer all these questions by 
investigating in more detail the newly developed thermo¬ 
dynamic theory for systems subject to a time-periodic 
perturbation [9|. The answers, which will be developed 
in further detail below, are the following. 1) As was al¬ 
ready anticipated by another thermodynamic argument 
in [l7j| . an Onsager matrix may be asymmetric but has 
to reduce to a symmetric form in the limit of zero dis¬ 
sipation. A crucial insight is that the reduced Onsager 
coefficients depend on the type of driving. The coeffi¬ 
cients change as one adapts the driving to move closer to 
the zero dissipation regime. It is in this very limit that 
the 2x2 Onsager matrix has to be symmetric. 2) From 
the expression of the entropy production in terms of all 
the fluxes in response to every possible Fourier mode, we 
find that zero dissipation implies that all fluxes vanish, 
and hence also all reduced or macroscopic fluxes. In par¬ 
ticular, any reduced 2x2 Onsager matrix will become 


symmetric in this limit. Hence this property is a result of 
the second law, provided its full impact for time-periodic 
perturbations is assessed. 3) The general implication on 
any Onsager matrix is that its determinant has to go to 
zero in the zero dissipation limit. For matrices of order 
higher than two, it does not need to reduce to a symmet¬ 
ric form. 

To proceed to these answers, we turn to the stochastic 
thermodynamics 18-2^ for a system described by a set 
of discrete energy levels subject to time-periodic mod¬ 
ulation. At this point, we do not need to specify the 
origin of this modulation. Suffice to say that the mod¬ 
ulation could be produced by several sources, which are 
supposed to be non-dissipative and thus describe the ex¬ 
change of work between the system and these sources. 
The system is furthermore in contact with a heat bath 
at temperature T, which can induce transitions between 
the different energy levels, entailing the exchange heat 
with the reservoir. We suppose, without loss of gener¬ 
ality, that this temperature is not modulated, but our 
conclusions apply equally well to, for example, thermal 
machines. We denote by tk the energy of level k. Its 
perturbation by a general time-periodic signal of period 
T is represented in terms of its Fourier decomposition: 


efc(t) — Cfc + EE {k,n,a)9{n^a) 1 


( 8 ) 


n—l (T—s,c 


with 


9{n,s){t) = sin 



9{n,c){t) = cos 



• (9) 


Ae(k,n,< 7 ) is the amplitude of the perturbation, applied 
to energy level fc, with Fourier mode n, and the a index 
referring to whether it is a sine (s) or cosine (c) perturba¬ 
tion. The rate of entropy production, averaged over one 
period (and still denoted, by slight abuse of notation, as 
S), is given by: 




/ _ 1 r - K 

Jo T TT ‘iJo T ‘IT’ 

( 10 ) 

where Q and W are the rate of heat and work to the 
system. In writing Eq. m, we assume that the system 
is in a periodic steady state. Its energy U returns to 
the same value after each period, which, combined with 
the first law, gives AU = Q + W = 0 for the integral 
over one period. The state of the system is described 
by a probability distribution p(t) = {pfe(t)}, with Pk(t) 
the probability to be in energy state k at time t. This 
distribution obeys the following master equation: 


p{t) =W{t)p{t). (11) 


Work on the system corresponds to energy dispensed 
upon moving an occupied energy level. In the present 
stochastic content, the rate of work is thus given by 
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W(t) = ek{t)pk{t). Combined with Eqs. ([8]) and (ITOt . 
this leads to the “familiar” expression for the entropy pro¬ 
duction as a sum of forces times fluxes Jq, (using the 
compact notation a = {k,n,a)): 

5 = ( 12 ) 

a. 

Fa = ^^, Ja = ^j dtg^^ri,a){i)Pk{t). (13) 

In the regime of linear response, the fluxes are linear 

functions of the forces: 

Ja = ^ LapFp, (14) 

/3 


with the Onsager coefficients given by: 




dJa 

dFp 


F=0 


(15) 


with a = (k, n, a) and /3 = {I, m, p). Here k and I refer to 
energy levels, n and m to Fourier modes, and a and p to 
the choice of a sine or cosine perturbation. Via a short 
calculation, cf. and supplemental material, one finds 
that the Onsager coefficients are the sum of an adiabatic 
and non-adiabatic contribution: 


The contributions per frequency, L„, can further be split 

into two contributions L„ = FlJn \ where is di- 

( 2 ) 

agonal in a and p, while L„ is anti-symmetric and hence 
does not contribute to the entropy production. Explic¬ 
itly: 







1 0 


0 1 
-1 0 


(T,p 


( 20 ) 


K%,i) = 




Mil n = 4 


eq 

p^Tin 






( 21 ) 


We are now ready to discuss the implications of zero 
entropy production. It requires thermodynamic forces F 
obeying FLF = 0. From the explicit structure of the 
Onsager matrix L, one finds: 


= ( 16 ) 
p^ipl^-Sk,) /•%,. ... 

F{k,m,a),(Un,p) = - ^(OS'(n.p) (O> 

eq !•% 

Triad _ t^l I 7, 

-^{k,m,<T),{l,n,p) ^ y 

Ik exp (1; - p'=‘^). (17) 

1; is the vector (0,..., 1,..., 0) with the value 1 on the I th 
position, = W(t)|p^ _Q is the time-independent un¬ 

perturbed evolution operator and the corresponding 
equilibrium distribution: 


/»LXJ 

/ dT - t) 

'n 


eq 

Pk = 



(18) 


'^FaLapFp ={) F(^k,n,<y) = F(l,n,a) V(fc,/)• (22) 

This can be seen by inspection from Eq. (EID by decom¬ 
posing the part of F, related to the energy levels, into 
the left eigenvectors of Assuming that this ma¬ 

trix is irreducible, the unit vector 1 = (1,!,...,!) is the 
only left eigenvector with a non-negative (zero) eigen¬ 
value (see supplemental material for more details). We 
conclude that zero entropy production is only compati¬ 
ble with a global modulation of the energy levels. Note 
that as a result the probability distribution p has to be 
time-independent. Hence, we expect that all fluxes will 
also be zero. Indeed, one immediately verifies that ther¬ 
modynamic forces of the above type are zero eigenvectors 
of the Onsager matrix: 


The above integrals can be performed by inserting the 
expressions given in Eq. ®. The adiabatic contribution 
describes the regime of an infinitely slow perturbation. 
It is anti-symmetric and hence does not contribute to 
the entropy production. Concentrating further on the 
non-adiabatic contribution, for which we use the matrix 
notation L, we first observe the “Curie principle” for 
time periodic variation: contributions from different fre¬ 
quencies do not mix. The Onsager matrix L can thus 
be written as a block diagonal matrix or as a direct sum 
over all frequencies n: 

L = ©„L„. (19) 


F(k,n,<7) = F(pri,a) 'd{k,l) ^ Ja = ^ Ta/jF/J = 0. (23) 

p 

This follows again from Eq. (EU, with the observation 
that 1 is also the right eigenvector of both and 

with eigenvalue zero. 

We finally return to the thermodynamic puzzle. As an 
illustrative example, we investigate on the basis of the 
preceding analysis the transformation between different 
sources of work i. The thermodynamic force F^, describ¬ 
ing the intensity of modulation of a given energy level, 
is due to the compound effect of all work mechanisms, 
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which we denote by the index i: 

i^a=^Fz,a. (24) 

i 


Introducing the thermodynamic forces Fi , prescribing the 
overall intensity of the driving i, and the corresponding 
reduced (coarse grained or macroscopic) fluxes: 


J^ = 


J2a Pi,aJa 


(25) 


one readily sees that the entropy production m assumes 
the reduced form S = The corresponding re¬ 

duced Onsager coefficients are obtained by the combina¬ 
tion of Eqs. (ITU) and (1^ : 

= E ( 26 ) 

a,/9 ^ 

We now make the crucial observation that, in contrast to 
the detailed Onsager coefficients Lq,/?, the reduced co¬ 
efficients Lij do depend on the thermodynamic forces 
Fa- Hence, zero dissipation, which imposes as discussed 
above conditions on these forces, will have an impact 
on the properties of the reduced Onsager matrix. In¬ 
deed, in the limit of zero dissipation, one concludes from 
Eq. (1251) that the reduced fluxes Ji go to zero since this 
is the case for the Jq,. Referring to our earlier discus¬ 
sion for a 2 X 2 Onsager matrix, we conclude that the 
latter has to become symmetric in the zero dissipation 
limit. An explicit illustration is given in Fig. [T] for a 
two-level system modulated by two work sources. More 
generally, the determinant of the Onsager matrix has to 
vanish (cf. supplemental material for an example of a re¬ 
duced 3x3 Onsager matrix, which is clearly asymmetric 
but has zero determinant). These properties can also be 
verified directly from the explicit expression Eq. dUl). 

In conclusion, we showed that zero dissipation implies 
the vanishing of all fluxes, implying in turn that any 
Onsager matrix must have a zero determinant. In partic¬ 
ular, a 2 X 2 Onsager matrix has to become symmetric. 
It should pose no problem to verify this prediction in 
a linear response experiment, involving time-periodic 
perturbation or magnetic fields, by moving closer to the 
zero dissipation regime. While the analysis given here 
does not readily apply to magnetic fields, we surmise 
that an analogous explanation will hold, probably 
requiring a more detailed thermodynamic analysis of 
electromagnetic phenomena. Finally, we mentioned that 
our analysis applies equally well to thermal machines. 
In fact, the illustrative example from provides a test 
case for our predictions: one easily verifies that heat and 
work flux become proportional to each other in the limit 
of reversible, i.e. Carnot efficiency, while the Onsager 
matrix becomes symmetric, cf. Eqs. (74) and (72) in this 
paper [^ . 



Figure 1. A two level system is modulated by two work 
sources 1 and 2, each operating on a different energy level 
of the system, namely ei(t) = ei -|- sin(27rt) -|- cos(27rt) and 
£2(1) = £2 -|- sin(27rt) -|-7cos(27rt), respectively. The param¬ 
eters are chosen such that = 3/4, = 1/4, T = 1, 

and nonzero eigenvalue of equal to —1. Main figure: 

Ji (blue) and J2 (green) are the powers of engines 1 and 
2, respectively, S (red) the entropy production, and lOdetL 
(orange) the determinant of the Onsager matrix. 7 > 1 and 
0.16 < 7 < 1 correspond to different operational regimes with 
2, respectively 1, functioning as driving and 1, respectively 2, 
as load. The corresponding efficiencies are 7 = —fi Ji/(f2 J2) 
and rj = —F2J2/{F\J\). For 7 < 0.16 both fluxes are larger 
than zero, and the efficiency is given by 7 = 0. Upper and 
lower inset: Onsager asymmetry Li2/i2i and efficiency 7 in 
function of the parameter 7. Zero entropy production with 
reversible efficiency 7 = 1 is reached for 7 = 1, in which case 
the Onsager matrix becomes symmetric. 
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I. ADIABATIC AND NON-ADIABATIC ONSAGER COEFFICIENTS: SUPPLE¬ 
MENTAL MATERIALS 


The Onsager coefficients are given by: 




dF(i^n,p) 
1 


F=0 




dpk{t) 


dF, 


{l,n,p) 


F=0 


The master equation; 


p(t) = lV(t)p(t), 

can be expanded in the detailed thermodynamic forces: 


with 


and 


Pjt) = JrP(i) 


d 


F=0 


= ^Wit) 


( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 


( 5 ) 


F=0 


The linear Eq. (|3]) can be solved explicitly. The contribution of the initial value disappears 
when reaching the time-periodic steady state: 


Pai.t) = 


dr exp - t)p^<^. 


( 6 ) 


It is convenient to introduce the instantaneous steady state of the full operator, which 

is reached when the system is driven adiabatically slow: 


W{t)p^\t) = 0, 


explicitly: 


Pfit) = 


-IkPl 

e 

H(i) ■ 


( 7 ) 


( 8 ) 


By expanding Eq. ([3]) in the thermodynamic forces, we arrive at the following relation: 


|y(0)pad(^) ^ = 0. 


( 9 ) 
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This allows to rewrite the quantities p^it) in terms of the unperturbed operator and the 
known adiabatic probabilities, Eq. ([8]): 


pjt) = - j <iT exp (t - t) 

= Pt(t) - r dr exp - r). 


( 10 ) 


where we integrated by parts with respect to r to obtain the second line. Returning to 
equation ([1]), we see that: 

L(k,m,a),{l,n,p) = j dt g(^rn,a){t)pl%t) - - j dt J dr lkg(m,a){t) tJ pl^{t - t) 


' ad 


= L 

where p'^{t) is given by: 

and Iz = (0,0,..., 1, ...,0). 


" nad 


I 


Ptn,p)it) = Pr9in,p)it) (p"'' - Iz) , 


( 11 ) 


( 12 ) 


II. EXPLICIT EXPRESSION OF THE ONSAGER COEFFICIENTS 


The adiabatic part of the Onsager matrix follows from: 
/•‘I 


'0 


di g(n,p)(j')9{m,a)(d') (1 ^cr,p)^n,m( 1) ’ ^Tn., 


namely: 


' ad 


^(/c,m,(j),(/,n,p) 


pT{pV - 4,z)(l - h^,p)h„,m(-l)‘^"’=7rn 


T 


For the non-adiabatic part, we hrst perform the integral over t\ 

I dt g(^jji^g'^{t') g(^n,s)ii "^) / dt g(^pi2^c)it') 9{n,c)^t t') 


= 26„ 


TT^n^ / 27 ™t\ 


T 


cos 


[—) 


(13) 


(14) 


(15) 


dt g(m,,s) (t')9{n,c) (t 


T = 


= 26, 


n,m 


dt g(^m,c)(j'')9(n,s)(d 

TT^n^ f 2'KnT\ 


(16) 
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To perform the r-intergral, we make a decomposition in right eigenvectors of 


P'»-1, = 5]o.A'. (17) 

i 

with A[ the right eigenvector corresponding to eigenvalue Aj. The corresponding left eigen¬ 
vector shall be denoted by A^. The eigenvectors are chosen to be orthonormal; A^A^ = 5ij. 
Note that the left eigenvector, Aq, corresponding to the supposedly unique zero eigenvalue 
of fulhlls Ag = 1 for all k, due to the orthonormality, and therefore: 


A!,p'’-A!,l, = 5;pr-1 = 0, (18) 

k 

which leads to ao = 0. Hence there is no contribution form the zero eigenvector. From 

exp = ^e^*"A[A'. (19) 


and the observation that all the other eigenvalues A* of are strictly negative, the 

integrals over r can be performed: 


dre^*'^ sin 


dre^*'^ cos 


/ 27rnrA 2‘l7rn 

V T ) ^ 47r2n2 + A^T^ ’ 

/ 27rnrA T^Aj 

V T ) ^ “47r2n2 +A2‘I2 


( 20 ) 


( 21 ) 


This leads to the following explicit expression for the non-adiabatic entropy Onsager coeffi¬ 
cients: 

2TTn ^ TA,- 


r nad 


^2 2 eq 

= ( (1 - 


47r2n2-|-A,?T2 °^’^47r2n2-1-A,?T2 / 

( 22 ) 

Now that all integrals have been executed, we can rewrite the result in matrix notation, by 
noting that : 

/(A,)A' = /(IV™)A). (23) 


III. ZERO-ENTROPY PRODUCTION 

To identify the conditions for zero entropy production, we decompose the thermodynamic 
force vector into a basis, consisting of vectors of the following form: 

= Fn®X\®(T, (24) 
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with Fn = In the ’natural’ basis in the frequency space, cr the ’natural’ basis in the 
sine/cosine space, and A[ the left eigenbasis of In this way, a general thermody¬ 

namic force vector, F can be written as: 

F' ^ ^(i,n,g-)-^(i,n,cr) • (25) 

i^n^cr 

We thus find for the rate of entropy production: 

S = FLF 

i^n^cr 

= _ ^ ^ 
i,n,a,k 

> 0. (26) 

This formula immediately shows that S' = 0 if and only if: 

b(^i^n,a) Si,ob{n,a): (27) 

i.e. the dependency on the energy level of the thermodynamic force is given by Ag = 1, 

hence the driving of the energy levels is uniform. Furthermore, it is clear that Ag is an 
eigenvector with eigenvalue zero of both and for all n. From this, Eq. fl27)) is a 
sufficient condition to have: 

J = LF = 0, (28) 

which can also be directly seen from Eqs. o and m- We conclude that zero entropy 
production implies zero detailed thermodynamic fluxes. 


n^\i 




[uKYpV 


IV. EXAMPLE: TWO LEVEL SYSTEM WITH THREE THERMODYNAMIC 
FORCES 


In this section we shall consider an extension of the example in the main text. Here we 
consider a three-level system, where the energy levels are driven as: 


27rf 


ei(f) = ei + sin ( — j + F(i,c) cos — J , 


2'Kt\ 


, , ^ f 2nt\ f 2'Kt\ 

ityt) = £2 + f( 2 ,.) Sin I — I + F( 2 ,„) ™S I ^ I . 


(29) 

(30) 
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= es + ^( 3 , 8 ) sin 


(31) 


In general, can be written as: 



+ F(3^c) cos 



|y(o) 


- k2pl'^ hpf hpT 

hpl^ -hpf - hpl^ k,pl^ 

k2p7 k^pl^ -k2pf - kspl'^ 


(32) 


We assume that there are three thermodynamic forces, which can be decomposed as: 


-^1,(1,s) — 2, -P'l,(l,c) — 1, -^ 2 ,( 2 ,s) — 1, -^ 2 ,( 2 ,c) — 7l) ^3,{2,s) — 72, -^3,(3,s) — 2, -^3,(3,c) — 72 

(33) 

where 71 and 72 are the variable parameters. The macroscopic forces are then the total 
amplitudes of the detailed forces: 

Fi = V5, Fa = + 7i, F 3 = 72 + y^4 + 7 |, (34) 

and the total detailed forces on the energy levels are found by summing up the contributions 
of the different systems: 

^il,s) = 2, -^(l,c) = 1, -^( 2 ,s) = 1+72 -^( 2 ,c) = 7l ^i3,s) = 2 +( 2,7 = 72- (35) 

The results are shown in Fig. [1] Note that the Onsager matrix is no longer symmetric in 
the zero dissipation limit, but the determinant still vanishes. 
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FIG. 1: Macroscopic fluxes Ji (blue), J 2 (green), J 3 (purple), the entropy production S (red) and 
10® det L (orange) in function of 7 with 71 = 72 = 7 - Parameter values: T = 1, = 0.7, p ^2 — 0-2, 

p^^ = 0.1, r = 1 A;i = 1, /c 2 = 0, ^3 = 2. Inset: L 12 /F 21 in function of 7 . 
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